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ABSTRACT 

We perform fully self-consistent stellar dynamical simulations of the accretion of 
a companion ("satellite") galaxy by a large disk galaxy to investigate the interaction 
between the disk, halo, and satellite components of the system during a merger. Our 
fiducial encounter begins with a satellite in a prograde, circular orbit inclined thirty 
degrees with respect to the disk plane at a galactocentric distance of six disk scalelengths. 
The satellite's mass is 10% of the disk's mass and its half-mass radius is about 1.3 kpc. 
The system is modelled with 500 000 particles, sufficient to mitigate numerical relaxation 
noise over the merging time. The satellite sinks in only ~ 1 Gyr and a core containing ~ 
45% of its initial mass reaches the centre of the disk. With so much of the satellite's mass 
remaining intact, the disk sustains significant damage as the satellite passes through. At 
the solar circle we find that the disk thickens ~ 60%, the velocity dispersions increase 
by Act ~ (10,8,8) km/s to (ctr, ct^, a z ) ~ (48,42,38) km/s, and the asymmetric drift 
is unchanged at ~ 18 km/s. Although the disk is not destroyed by these events (hence 
"minor" mergers), its final state resembles a disk galaxy of earlier Hubble type than 
its initial state, thicker and hotter, with the satellite's core enhancing the bulge. Thus 
minor mergers continue to be a promising mechanism for driving galaxy evolution. 

Subject headings: galaxies: evolution — galaxies: interactions — galaxies: structure - 
Galaxy: evolution — Galaxy: kinematics & dynamics — Galaxy: solar neighbourhood 



Hubble Fellow 

present address: Johns Hopkins University, Department of Physics and Astronomy, Baltimore, MD 21218 
Alfred P. Sloan Foundation Fellow, Presidential Faculty Fellow 



1. Introduction 

The business of simulating interactions between 
galaxies originated with Holmberg (1941!) who found 
that a close encounter could raise impressive and ob- 
servable tidal distortions and even thcrmalize enough 
orbital energy to result in capture. Although he was 
studying the clustering of galaxies, this work came 
long before the recognition of large-scale structure 
and it was revisited only rarely (e.g., by Pfleiderer 
& Siedentopf (1961); Pfleiderer (1963)) because of 
the cherished view of galaxies as "island universes" 
which rarely interact. However, a number of observed 
galaxies show signs of interaction (Arp 1966) and the 
subject was reintroduced around 1970, particularly 
by Toomre & Toomre (1972), who successfully mod- 
elled several peculiar galaxies with collisions. The 
expected frequency of encounters grew with the be- 
lief that galaxies are embedded in extensive, massive, 
dark halos which increase cross-sections for major col- 
lisions and enable the orbits of satellite galaxies to 
decay by dynamical friction (Tremaine 1981). Now, 
more than fifty years after Holmberg's illuminating 
efforts, the prevalence of the hierarchical clustering 
picture for the formation and evolution of structure 
makes the study of galaxy interactions a very active 
field and a wide range of phenomena are thought to be 
linked to one variety of galaxy interaction or another. 

Major mergers between spirals of comparable mass 
have received the most attention. The intermediate 
stages of such collisions explain some of the most 
spectacular objects observed, reproducing their messy 
profiles and extended tidal tails. Because the stellar 
component of the remnant relaxes toward a de Vau- 
couleurs profile, such mergers appear to drive evo- 
lution along the Hubble sequence toward elliptical 
galaxies (Negroponte & White 1983; Barnes 1988, 
1992; Hernquist 1992, 1993a). They may also trigger 
starbursts or nuclear activity when gas is driven to 
the remnant's centre as the progenitors coalesce (Mi- 
hos & Hernquist 1994a; Barnes & Hernquist 1995). 

Less spectacular but more common, mergers in 
which a disk galaxy accretes a smaller object ("mi- 
nor mergers") have analogous effects, namely, the 
stirring-up of disks and the generation of peculiar fea- 
tures (Quinn & Goodman 1986; Quinn, Hernquist & 
Fullagar 1993), bulge-building (a smaller step along 
the Hubble sequence toward Sa/SO (Schwcizer 1990)), 
and enhanced star formation or nuclear activity (Hern- 
quist 1989; Mihos & Hernquist 1994b). Specifically, 



the early stages of such interactions (or even grazing 
encounters) can induce spiral arms, bars, warps, or 
"bridges" and full merging can thicken and dynam- 
ically heat disks. The tidal stripping of a satellite 
can produce features such as long tails or counter- 
streaming groups and its accretion introduces a new 
stellar population into the disk. A galactic bulge may 
be formed or enlarged by any combination of (a) the 
satellite's core surviving to the centre of the disk, (b) 
gas being driven to the centre and converted to stars, 
or (c) disk stars being stirred up out of the disk plane. 
Combine this bulge-building with the smoothing of 
disk features which results from the increased veloc- 
ity dispersion and the primary galaxy looks more like 
an SO than it did before. For a summary of these 
ideas see the review by Barnes & Hernquist (1992) or 
any one of a number of conference proceedings (that 
of Wielen (1990) is comprehensive). 

Those classic, messy objects which spark observa- 
tional interest in major mergers have been known for 
decades. Interest in minor mergers has mostly devel- 
oped more recently because the subtler features at- 
tract less attention and can be much harder to de- 
tect and study. However, considerable observational 
data has accumulated. The Milky Way has a num- 
ber of satellite companions, as do most large galax- 
ies. There is evidence that some quasar activity could 
be associated with minor mergers (Bahcall, Kirhakos, 
& Schneider 1995). Grand design spiral structure 
is sometimes associated with small, close compan- 
ions: M51 is a classic example. A large fraction 
of disk galaxies are warped (Binney 1992). "Thick 
disks", hot, flattened components with kiloparsec 
scaleheights, have been discovered in our Galaxy 
(Gilmore & Reid 1983; Gilmore, Wyse, & Kuijken 1989) 
and others (Burstein 1979; van der Kruit & Searle 1981). 
The Magellanic Clouds and Stream have been known 
for a long time, but the discovery of the close, tidally 
stretched Sagittarius dwarf galaxy (Ibata, Gilmore, & 
Irwin 1994) and the recognition of a distinct, possi- 
bly accreted population of blue, metal-poor stars in 
the solar neighbourhood (Preston, Beers, & Sheet- 
man 1994) have generated recent excitement. Most 
recently, we have found that the bulge of Hickson 87a 
very closely matches the bulge produced in one of our 
simulations (Mihos et al. 1995). The parallel between 
the observational results listed here and the minor 
merger effects listed in the previous paragraph drives 
the current interest in discrete accretion events. 

Tremaine (1981) used the mean radial distribution 
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of satellite companions of large galaxies and sinking 
rates derived from dynamical friction in an isothermal 
halo to estimate the total mass accreted by a large 
galaxy over a Hubble time. He found that a typical 
large galaxy can consume 20%-30% of its own mass, 
very roughly, motivating several attempts to quantify 
the resilience of cold, thin disks to minor mergers: 
What are the implications of the chunk-wise accre- 
tion of so much mass? Quinn & Goodman (1986) 
measured the effect of the disk on these sinking rates, 
studying the dynamics of the disk-satellite interac- 
tion with analytic and restricted iV-body techniques. 
They noted strong heating of their simulated disks. 
Quinn, Hcrnquist, & Fullagar (1993) used A-body 
experiments to survey orbital parameters, focussing 
on such observable disk properties as density struc- 
ture and kinematics. They found that single accre- 
tions of about 10% of the disk's mass could, in the 
solar neighbourhood, double the disk's thickness and 
expand the velocity ellipsoid by fifty per cent (more 
in the radial direction), producing something similar 
to the Milky Way's thick disk. The analytic work of 
Toth & Ostriker (1992) also found disks to be quite 
fragile. They turned Tremaine's question around by 
using the observed abundance of undisturbed, cold, 
thin disk galaxies to limit the mean total mass ac- 
creted and, by inference, Slo ■ They concluded that 
the Milky Way cannot have accreted more than 4% 
of the mass inside the solar circle in the last 5 Gyr. 
Since f2 = 1 models predict that accretion continues 
at late times, they argued that < 1. 

Successes notwithstanding, these studies had the 
drawback of not being fully self-consistent, modelling 
the dark matter halo and sometimes also the satellite 
as rigid bodies. Modelling a component with inter- 
acting particles (thus making it responsive, or "live" ) 
increases both computational expense and noise and 
might seem frivolous for the halo which is, after all, 
invisible. However, rigid representations lack inter- 
nal degrees of freedom and thus the disk was the 
only sink for the satellite's energy in earlier inves- 
tigations. These studies were in some sense measures 
of the maximum damage that could be done to the 
disk, not of the typical response. Also, a live halo is 
able to propagate disturbances into the disk (Wein- 
berg 1995a), a pathway which is cut off when a rigid 
halo is used. 

Confidence in the results of Quinn et al. (1993) 
is undermined by several other technical limitations 
and compromises, apart from the use of a rigid halo. 



Their simulations used 32 768 disk particles and 4096 
satellite particles and were thus fairly noisy. Their 
satellite was initialized with some fraction of its mass 
not bound because the tidal radius was not taken into 
account. This mass was left behind at large radii and 
did not participate in the stirring of the disk. Also, 
their disk was abruptly truncated at the initial ra- 
dius of the satellite's orbit, interfering with the disk- 
satellite coupling. This allowed the satellite to or- 
bit several times at large radii before sinking, making 
sinking times unreliable but also allowing the satellite 
to shed even more mass at large radii, further reduc- 
ing its potency for scattering disk particles. 

Advances in computer technology now permit the 
use of an order of magnitude more particles than 
used by Quinn et al. (1993). Here we attempt to 
remedy the shortcomings of earlier studies with fully 
self-consistent N-body simulations large enough to 
eliminate numerical relaxation noise over the merg- 
ing time while including the dynamical interplay of 
all the components of the system. 

2. Models and Methods 
2.1. Galaxy Models 

Our models consist of exponential disks, trun- 
cated isothermal halos, and Hernquist model satel- 
lites (Hernquist 1990a) in the mass ratio 10:58:1. The 
disk and halo each contain 45% of the simulation par- 
ticles while the satellite contains the remaining 10%, 
an arrangement which provides better sampling of the 
luminous components. The specific profiles used are: 

Pdisk(^) = P(0)e- fl /' l cosh- 2 (z/zo), (1) 

e -r 2 /r 2 c 

PhaloW = K°) 1 + r2/7 2 ' ( 2 ) 

, . Ma 
P ^ {r) = 2nr{r + af (3) 

The satellite is truncated at its initial tidal radius 
with respect to the primary. Parameter values are 
listed in Table 0. 

Velocities are initialized from moments of the colli- 
sionless Boltzmann equation in a procedure described 
by Hernquist (1993b). For each dark matter particle, 
the velocity ellipsoid suitable for its location is calcu- 
lated from the moment equations and then its veloc- 
ity components are randomly selected from gaussians 
corresponding to that ellipsoid. For disk particles, 
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velocities are constrained to make the radial disper- 
sion proportional to the square root of surface density 
(<jr oc e r R l' lh \ Freeman 1993) and are normalized 
to a chosen value of Q(R Q ). The vertical dispersion 
is given by o\ — 7rGS(i?)z (the so-called isother- 
mal sheet). The azimuthal dispersion is given by the 
epicycle approximation (er^ = aj^K 2 /(Afl 2 )) and the 
azimuthal streaming velocity by the cylindrical mo- 
ment equations. The initial disk structure is shown 
in Figure 

This approach is not without its drawbacks since 
it only approximates a distribution function and thus 
does not initialize the particles in a true equilibrium 
configuration. When this initial state is allowed to 
evolve in isolation (i.e., without a satellite) it rapidly 
shifts to an equilibrium configuration whose subse- 
quent evolution is governed only by numerical relax- 
ation (simulation noise). However, this shift can be 
large. Figure || illustrates this by superposing the ini- 
tial disk structure on the structure after several time 
steps. Note in particular that the central dispersions 
have dropped by about twenty per cent. 

For our purposes this technique is adequate be- 
cause the galaxy is about to be stirred up by an in- 
falling neighbour. As long as we are careful to mea- 
sure "heating" and "thickening" relative to the coeval 
isolated disk galaxy rather than the initial state, our 
analysis will be sound. On the other hand, studies 
involving isolated disk galaxies (e.g., instability stud- 
ies) would require better initial conditions, both for 
tighter control of initial structure and less disturbance 
from transients. Improvements can be obtained from 
higher order moments but for future work the semi- 
analytic distribution functions of Kuijken & Dubinski 
(1995) may offer a more direct approach. 

2.2. Encounter 

We choose to focus on an encounter in which the 
satellite starts on a circular, prograde orbit with a 
radius of six disk scalelengths inclined 30° with re- 
spect to the disk plane. This encounter was found by 
Quinn et al. (1993) to illustrate most of the impor- 
tant phenomena. Note that the satellite is not grown 
adiabatically in this orbit; it is simply switched on at 
t = 0. Although the initialization procedure tries to 
account for the effect on the satellite by truncating it 
at its tidal radius, the disk is thrown somewhat out of 
equilibrium by the sudden change in potential at its 
edge. However, as described above and demonstrated 
in Figure 0, the disk is initially out of equilibrium any- 



way, and the transient behaviour within 7? = 15 kpc is 
identical to that for the isolated galaxy shown in Fig- 
ure ^. Nevertheless, as computers become ever faster, 
it will be prudent to invest some cpu time in moving 
the starting radius out farther for the sake of more 
realistic orbital evolution, tidal stripping, flaring and 
warping of the disk, et cetera. 

For the sake of discussion, units are scaled to 
the Milky Way such that the disk has scalelength 
h = 3.5 kpc and mass M d = 5.6 x 10 10 Af Q (Bah- 
call, Schmidt, & Soneira 1983). Other parameters are 
given in Table [l]. The evolution was followed with a 
treecode (Barnes & Hut 1986; Hernquist 1987, 1990b) 
using terms up to quadrupole order, a tolerance pa- 
rameter ("opening angle") 9 = 0.7, and a timestep 
At ~ 2 Myr. Energy was conserved to 0.1% over 
1.2 Gyr and to 0.2% over 2.5 Gyr. Simulations which 
include hydrodynamics have also been studied and are 
reported separately (Hernquist & Mihos 1995). 



Quantity 


Value 


Disk 




Mass M d 


5.6 x 10 10 M Q 


Scalelength h 


3.5 kpc 


Scaleheight zo 


700 pc 


Softening ea 


140 pc 


Toomre Q(R Q ) 


1.5 


N d 


0.45A a 


Halo 




Mass M h 


3.25 x 10 11 M Q 


Scalelength 7 


3.5 kpc 


Cutoff Scale r c 


35 kpc 


Softening eh 


700 pc 


N h 


0.45A^ a 


Satellite 




Mass M s 


5.6 x 10 9 M Q 


Characteristic Radius a° 


525 pc 


Softening e s 


70 pc 


N s 


0.1AT a 


Orbit i? init 


21 kpc 



Table 1: Simulation Parameters 



a A r = 500 000 particles for the largest simulation. 
b The half-mass radius is (1 + v2)o. 

2.3. Quality Control 

Modelling a halo (or any other component) with 
particles does have a disadvantage in that the po- 
tential is less smooth than a rigid, analytic halo's 
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and scattering is thus enhanced. Compared with real 
galaxies, simulated galaxies have a relatively small 
number of relatively massive particles whose wide, 
deep potential wells are able to deflect one another 
even if they pass at a significant distance. This excess 
scattering, called 2-body relaxation noise, is to some 
extent suppressed by "softening", artificially reduc- 
ing the force between close particles to prevent strong 
collisions. However, random clumping ("shot noise") 
creates potential fluctuations on all scales and soft- 
ening does not help with scattering off large clumps. 
The contrast between random clumps and the back- 
ground is reduced as more particles are used; if we 
are to detect small but legitimate structural and kine- 
matic effects in our simulations, we must reduce the 
excess scattering as much as possible. To quantify 
this, we ran several simulations of the disk galaxy 
in isolation (i.e., without the satellite). Figure || 
shows the evolution of the disk thickness and veloc- 
ity ellipsoid at the solar radius in simulations with 
N = 45 000, which is roughly the size of the Quinn 
et al. (1993) simulations, TV = 90 000, N = 225 000, 
and N = 450 000. The largest run experienced very 
little change over the sinking time (~ 1 Gyr) so we 
expect any evolution we observe in our satellite en- 
counters to be signal rather than noise. 

Shot noise has also been observed to have a global 
manifestation. Disk models with realistic profiles 
and velocity dispersions tend to be quite lively in 
that they amplify density perturbations which can 
then feed into instabilities and produce spiral or bar 
features (Toomre 1981; Binney & Tremaine 1987; 
Sellwood 1989; Hernquist 1993b; Weinberg 1995b). 
This is exacerbated in models with live halos because 
clumps of massive particles create wakes in the disk. 
Since the halo dominates the potential and is the most 
poorly-sampled component, it is the principal source 
of shot noise. This is apparent in Figure ||, which 
shows the growth of the bar (to = 2) mode in our iso- 
lated galaxies and in a run with 45 000 disk particles 
and a rigid (noiseless) halo whose bar mode seems 
not to grow at all. Because the contrast between the 
random clumps and the underlying smooth distribu- 
tion is smaller when N is larger, the instability is 
seeded at a lower amplitude and thus sets in at later 
times in larger simulations. Again, it is apparent that 
noise has been suppressed for the relevant timescale 
in the largest simulation and that any significant bar 
growth in our satellite encounter can be assumed to 
be induced by the satellite. 



Thus we use the 450 000-particle disk galaxy with 
a 50 000-particle satellite to model our chosen en- 
counter. Figure [5] and Figure ^ demonstrate the need 
for large N to beat down noise over long timescales: 
i-million particles is only adequate for about two bil- 
lion years. Many more halo particles must be used for 
examining phenomena with timescales longer than a 
few gigayears, such as retrograde encounters. 

3. Results 

The basic sequence of events is illustrated in Fig- 
ure [5] which shows both face-on and edge-on views 
of the disk and satellite at regular intervals of about 
125 Myr. The satellite loses most of its vertical mo- 
tion while completing only ~ 1.5 orbits, settling into 
an orbit coplanar with the disk. This orbit then de- 
cays quite rapidly. The satellite sheds mass all along 
its orbit but its core survives and arrives at the centre 
of the disk. 

In this paper we discuss the decay of the satellite 
orbit and the satellite's effect on the disk. However, 
the disturbance continues to evolve for some time af- 
ter the merger is complete. The evolution of the cen- 
tral regions at late times is presented in a separate 
paper (Mihos et al. 1995), but here we give a brief 
description. Figure || illustrates the disk's global re- 
sponse to the satellite. After the merger (Fig. ^ last 
frame), full axisymmetry is not restored: a bar has 
been induced and it drives further evolution in the 
disk's inner two scalelengths. Note that this bar is 
induced by the satellite, not by amplified noise. Be- 
cause the bar is unstable to vertical bending (Raha 
et al. 1991), it flexes and eventually kicks material 
out of the disk plane, generating a small, flattened 
bulge with an X- morphology (Mihos et al. 1995). For 
our purposes, however, it is sufficient to note that, 
apart from changes associated with the bar, the disk 
structure does not change much, postmerger. 

3.1. Orbital Decay 

Figure || shows various aspects of the orbital decay. 
The upper left panel shows cylindrical radius versus 
time. This quantity does not decrease monotonically 
because the satellite is initially orbiting in an inclined 
plane. Once the satellite has settled closer to the 
disk plane, there is a knee in the R versus t curve 
and the orbit decays very rapidly. The upper right 
panel shows the orbit from the North galactic pole 
and reveals that the satellite falls most of the way to 
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the centre in only two orbits (R ~ 5 kpc). The lower 
left panel shows the altitude versus time and the lower 
right shows altitude versus radius, demonstrating that 
the satellite settles most of the way to the disk plane 
while still at large radii. 

In an attempt to quantify the importance of the 
halo compared with the disk over the course of the 
encounter we calculate the torque on the satellite due 
to the disk and halo separately (Fig. [?], upper panel). 
The two torques are comparable initially but the disk 
dominates overall. A rough integration of the curves 
in the upper panel of Figure [?] shows that ~ 75% of 
the total time-integrated torque comes from the disk. 

This is attributable to a resonant enhancement of 
dynamical friction which occurs because a satellite on 
a prograde orbit moves with the nearby disk particles 
and thus interacts strongly with them. The interac- 
tion dynamics were discussed in detail by Quinn & 
Goodman (1986), but basically leading disk particles 
are pulled back by the satellite so that they lose an- 
gular momentum and drop into lower orbits, clearing 
a gap in the disk in front of the satellite (Fig. ||). 
Trailing particles gain angular momentum and move 
into higher, slower orbits and thus do not overtake 
the satellite. Once this configuration develops, with 
the density higher behind the satellite than before it, 
the disk torque increases dramatically and the satel- 
lite quickly sinks to the centre. The lower panel of 
Figure superposes the bar mode strength, the disk 
torque, and the orbital decay to show that these co- 
incide. Note that this mechanism transfers most of 
the satellite's orbital energy into the potential energy 
of the disk. If, when the satellite pulls a trailing disk 
particle forward, promoting it to a higher orbit, the 
particle were to move from one circular orbit to an- 
other and the rotation curve were exactly flat, then 
the particle's kinetic energy would not change a bit 
and all the energy flowing from satellite to disk would 
appear as potential energy. 

Clearly, the process is aided by the settling of the 
satellite toward the disk plane, although the rapid 
radial decay of the orbit actually begins before the 
vertical motion has damped out completely (Fig. ||). 
The process also relies on the approximate circular- 
ity of the orbit in order for the velocity of the satel- 
lite to match the velocities of the disk particles. The 
global nature of the response presumably requires an 
approximately flat rotation curve to make the reso- 
nant region extensive. The strength may depend on 
the ease of excitation of the bar mode; however, a 



run with a flattened Hernquist model bulge (semi- 
axes a — 0.7, c = 0.35 kpc) yielded almost identical 
sinking time and heating. 

The halo, whose particles do not collectively orbit 
with the satellite, is unable to undergo strong collec- 
tive interaction and so cannot contribute such a large 
torque. Of course, in real life the satellite does not 
start at the edge of the disk. The torque on a satellite 
at large distance will be dominated by the halo and 
the satellite will take a long time to sink toward the 
disk where the process described above can set in. In 
fact, the disk's contribution to the torque in the early 
stages of the decay of a prograde orbit can oppose the 
halo's (Quinn & Goodman 1986): the satellite can 
gain angular momentum from the disk in a manner 
analogous to the Moon's gain of angular momentum 
from the tidal bulge it raises on the Earth. This makes 
the truncation of the disk by Quinn et al. (1993) prob- 
lematic. Their Figure 11 shows orbital evolution sim- 
ilar to that seen here but delayed. The knee in their 
R versus t curve comes at ~ 2 Gyr and the total sink- 
ing time is ~ 2.5 Gyr. 

We test our simulations for convergence by run- 
ning four simulations of our encounter with 50 000, 
100 000, 250 000, and 500 000 particles, respectively. 
The upper panel of Figure ^ shows the orbital decay 
curves. They converge well and, since our simulations 
are fully self-consistent, we conclude that we finally 
have reliable estimates for the sinking time. It is not 
very long: at 1 Gyr, it is only ~ 10% the age of the 
Milky Way's disk. (For comparison, we overlay the re- 
sult from a simulation identical to our largest in every 
respect but for a rigid halo. The sinking time is some 
50%-80% longer (though still much shorter than in 
Quinn et al. (1993)) which testifies to the importance 
of a self-consistent treatment.) Even more dramatic is 
a coplanar (i = 0°), prograde encounter in which the 
satellite sinks from six scalelengths in 0.6 Gyr (Fig. pi 
lower panel) . Evidently, once a satellite arrives in the 
disk on a prograde orbit, it has only a short time 
to live. This is in stark contrast to retrograde and 
polar orbits. We tried evolving them for ~ 3| Gyr 
and what little was left of the satellite core had just 
reached the galactic centre (Fig. 3, lower panel). As 
is apparent from Figures || and i, noise becomes a 
problem on these timescales. 

3.2. Disk Structure 

Figure ^ shows the disk thickness and velocity dis- 
persions as functions of cylindrical radius. The solid 
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lines show the structure at ~ 1.2 Gyr (which al- 
lows ~ 200 Myr for things to settle after the satel- 
lite reaches the centre) and the dashed lines show the 
structure of the coeval isolated galaxy. The entire disk 
thickens with respect to the isolated galaxy but much 
more so at large radii where the satellite still had sig- 
nificant vertical motion: more than 200% thickening 
is seen beyond R ~ 15 kpc, as opposed to ~ 50% at 
R ~ 5 kpc and only ~ 10% at the centre. In contrast, 
the disk thickens much more uniformly in the copla- 
nar encounter but only by about 20%, or 100 parsecs. 
The velocity dispersions in Figure || increase at least 
10 km/s at all radii. At the centre, an increases by 
50 km/s and a z increases by 30 km/s, which is in- 
teresting because the central thickness is essentially 
unchanged. This reflects the increased depth of the 
central potential due to the satellite core (see § 3.3). 
Note that the solid lines in Figure || are more ex- 
tended than the dashed: the disk spreads radially as 
well as vertically, storing much of the energy it gains 
as potential energy. 

The evolution of these quantities at the solar cir- 
cle (8.0 kpc) is shown in Figure [Tc|. Recall that these 
quantities are azimuthally averaged and that the disk 
is not axisymmetric at intermediate times: the num- 
bers calculated for these times should not be trusted 
absolutely. However, trends in the disk's behaviour 
are more reliable. The figure shows that the satel- 
lite's effect is felt in the radial dispersion first and in 
the other quantities only when the satellite reaches 
the solar circle. This pattern occurs in all the simu- 
lations and so should be robust. The increase in the 
radial dispersion is due to azimuthal averaging when 
there is radial streaming in the bar. It is thus not 
surprising that ctr increases early because the bar de- 
velops before the satellite reaches the solar circle. The 
asymmetric drift also responds to the passage of the 
satellite but settles back to about 18 km/s so that, 
within the noise, it has hardly changed at all. The 
other quantities settle but not to their original val- 
ues. In the end, the disk thickens by about 60%. The 
radial dispersion goes from about 35 km/s to about 
48 km/s, the azimuthal from about 32 to about 42, 
and the vertical from 28 to 38. Compare these fi- 
nal values with those from the isolated disk: around 
t ~ 1.2 Gyr, it has actually shrunk about 3%, heated 
to (38, 34, 30) km/s, and has an asymmetric drift of 
16 km/s. The net heating attributable to the action 
of the satellite is thus (10, 8, 8) km/s. 

The satellite acts to heat and thicken the disk by 



scattering disk particles. Because the satellite is con- 
centrated and moving, disk particles on close trajecto- 
ries can encounter slightly different parts of the satel- 
lite potential and can be deflected in slightly different 
directions. This creates a greater variety of local tra- 
jectories, i.e., the velocity dispersion increases. So, 
while the satellite's energy is absorbed into the po- 
tential energy of the disk (and the stripped satellite 
material), the heating and thickening of the disk occur 
more by pure scattering, the conversion of disk orbital 
energy to disk thermal energy. Figure [l(] supports 
this view by indicating that heating and thickening 
of the disk at the solar radius occur as the satellite 
passes through that radius and are thus essentially 
local processes. 

3.3. Satellite Remnant 

That the satellite's core arrives as a distinct lump 
in the galactic centre is apparent from the integrated 
mass distribution and surface density of the satellite 
remnant (Fig. [ll]). About 45% of the satellite's mass 
lies within 2 kpc of the centre at the end, enhancing 
the peak surface density by a factor of about 2.5 with 
respect to the disk. Contrast this with the retrograde 
case in which the satellite barely survives. 

Estimating mass loss with a simple density crite- 
rion based on the "tidal radius" for a body in orbit 
about another (Binney & Tremaine 1987) is inade- 
quate for the process we are studying. Obviously, 
any treatment which approximates the disk as a point 
mass when the satellite is orbiting within the disk 
is inappropriate, but such estimates have been used 
(e.g., by Toth & Ostriker (1992)) to relate mass loss to 
orbital radius. To see the potential for error, consider 
that the density ratio criterion does not make use of 
the spin of either body, hence cannot distinguish be- 
tween prograde and retrograde events, and thus pre- 
dicts the same mass loss for both. Although our ret- 
rograde simulation is unreliable for its fine structure 
because of the growth of noise over the long sinking 
time, its bulk properties are sound. The initial con- 
figuration is identical to that of our main simulation 
except for the sign of the orbital angular momentum; 
the mass loss history is quite different. By the end of 
the simulation, the satellite core has been pared down 
to only 10% of the initial satellite mass (Fig. [ll], upper 
panel). If half the satellite survives to the centre in 
the prograde case while the satellite is almost entirely 
disrupted in the retrograde case then a mass loss es- 
timator which does not distinguish these cases is all 
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but useless. Since it takes into account neither the 
angular momentum of the objects nor the timescale 
of interaction, the density criterion does not capture 
the physics involved and is misapplied in this prob- 
lem. Given that a satellite's mass loss determines its 
potency for scattering disk particles, any study rely- 
ing on simple analytic estimates of the mass loss must 
be interpreted very carefully. 

As for the material which is stripped from the 
satellite, it is distributed in a thickened, flared disk, 
similar to the final distribution of disk material but 
thicker. Because of its low surface density, this mate- 
rial would not stand out in an external galaxy; only 
the core, visible as a central brightness enhancement 
or small bulge (the peak in the dotted curve in the sec- 
ond panel of Figure |ll|), would be conspicuous. We 
note here that Hickson 87a, a galaxy whose luminosity 
structure matches our model quite closely, has such a 
peak (Mihos et al. 1995). Figure [l2| shows the struc- 
ture of the satellite remnant and disk combined, along 
with the disk and satellite separately for comparison. 
Here again the core stands out as a distinct compo- 
nent. Otherwise, the satellite remnant is structurally 
and kinematically like a thick, flared, hot disk. There 
is little difference between the disk structure and the 
combined structure (except at very large radii) be- 
cause the surface density of the remnant is just too 
small. Near the solar circle, the satellite remnant is 
somewhat thicker and hotter than the disk material 
but it does not lag in its rotation (within the noise). 

Since the satellite material does not alter the com- 
bined structure much, especially at Rq, it is not 
clear that it could be distinguished without a spec- 
tral signature (e.g., the blue, metal-poor population 
discussed by Preston et al. (1994)). At intermedi- 
ate stages in an accretion event, however, a satellite 
would be more noticeable (e.g., the Sagittarius dwarf 
galaxy (Ibata et al. 1994)) and streams and moving 
groups can persist for more than a gigayear (John- 
ston, Spergel, & Hernquist 1995). Of course, the 
counter-streaming material left by a retrograde en- 
counter would be easily distinguished, as would some- 
thing resembling the counterrotating cores studied by 
Balcells & Quinn (1990). 

4. Discussion 

4.1. Is our puffed up disk a "thick disk"? 

The Milky Way's thick disk is observed to be about 
twice as thick as our model's at the solar circle. Mor- 



rison (1993) quotes for the (metal-strong) thick disk 
a scaleheight of about 1 kpc, a z ~ 40 km/s, and an 
asymmetric drift of about 30 km/s. The velocity el- 
lipsoid is observed to be (<7r, cr^,,cr z ) = (63 ± 7, 42 ± 
4,38 ± 4) km/s by Beers & Sommer-Larsen (1995). 
Our simulated disk has a smaller asymmetric drift al- 
though a more eccentric satellite orbit might induce 
more lag. Its dispersions are (48, 42, 38) km/s, plus or 
minus a few kilometres per second, a bit cool radially 
but otherwise in good agreement. It is important not 
to overinterpret this comparison, however. The ex- 
act values we obtain depend sensitively on the choice 
of scaling (Table @). Moreover, our simulation starts 
with a fully formed galaxy and thickens its entire disk, 
whereas the Milky Way has a low-mass thick disk with 
an embedded high-mass thin disk. Thus the thick 
disks in the model and Milky Way are not strictly dy- 
namically equivalent and so not directly comparable. 
The Milky Way's thick disk is, after all, ancient, so 
while simulations of this sort help us study particular 
processes with few competing effects, honest compari- 
son with observations ultimately requires that they be 
studied in the context of the formation and evolution 
of the galaxy as a whole. 

The discreteness of the Milky Way's thick disk 
with respect to the thin disk argues for at least 
one significant ancient merger event. A tantaliz- 
ing hint for another is the suggested substructure 
in the age-velocity dispersion relation discussed by 
Freeman (1993), analyzing a figure from Edvardsson 
et al. (1993). Freeman proposes that the relation con- 
tains these three domains: stars younger than 3 Gyr 
with <j z ~ 10 km/s, stars between 3 and 10 Gyr with 
<r z ~ 20 km/s, and stars older than 10 Gyr with 
a z ~ 40 km/s (the thick disk). Whether the tran- 
sition between the first two domains is really discrete 
is not clear but the figure (Figure 2 in Freeman (1993), 
Figure 16b in Edvardsson et al. (1993)) at least hints 
at the possibility of a second, smaller, accretion event 
about 3 Gyr ago. 

The thickness of our simulated disk increases at 
large radii because of the inclined satellite orbit but 
is essentially uniform near the solar circle. At late 
times, however, the bar's vertical instability causes 
additional thickening (to (\z\) ~ 1.1 kpc) at radii 
around 1 scalelength. Interestingly, an increase in the 
Milky Way's thick disk scaleheight inside the solar cir- 
cle is one interpretation Morrison (1993) discusses for 
her data, though it is not strongly supported. Any 
detection of scaleheight varying with radius is likely 
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to come from external galaxies, but observers face the 
obstacle (which we do not) of having to subtract thin 
disk and bulge profiles before searching for thick disks 
in noisy residuals. In most instances, the issue is 
whether or not a thick disk has been detected, not 
what its detailed structure is, so the data is generally 
fit with a single thick disk scaleheight (van Dokkum 
et al. 1994; Morrison, Boroson, & Harding 1994). The 
new generation of large telescopes and large CCDs 
is expected to bring the study of thick disk profiles 
within reach. 

Given the number of parameters available for twid- 
dling (e.g., satellite mass, density, orbital inclination, 
and orbital eccentricity) it seems plausible that we 
could arrange a suitable thick disk by satellite accre- 
tion. Of course, there is no thin disk left, a result 
which motivates the study of models which include 
gas (Mihos & Hernquist 1994b; Hernquist & Mihos 
1995). The satellite can puff up both the gas and the 
existing stellar component but the gas can radiate 
away its energy and resettle to form a new thin disk, 
leaving a thick disk composed of a necessarily older 
stellar population. The degree to which the thick and 
thin disks would be distinct depends on the timescale 
for resettling and the ability of the gas to form stars 
as it resettles. 

4.2. Is satellite accretion a good way to form 
bulges? 

Certainly small nuclear concentrations can be pro- 
duced if the satellite core reaches the centre. As 
pointed out in § 3.3, Hickson 87a sports such a fea- 



ture. There will also be a contribution from that part 
of the gas which is driven to the centre (Mihos & 
Hernquist 1994b; Hernquist & Mihos 1995). It is not 
so clear that anything extended could be produced. 
A more extended satellite would simply be stripped. 
A much more massive or dense satellite would be too 
destructive for the disk. The bar's vertical instability 
kicks disk material up out of the disk plane to pro- 
duce a modest extended bulge (Mihos et al. 1995), 
though depending on viewing angle it may be a pecu- 
liar one (boxy, peanut-, or X-shaped). That our ini- 
tially bulgeless simulation ends up looking very much 
like Hickson 87a (an SO pec) is encouraging but gas in- 
flow may interfere with the bar process because steep- 
ening of the rotation curve tends to dissipate bars 
(Hasan & Norman 1990). Even without gas inflow, 
a more bar-stable model might have quickly returned 
to axisymmetry instead of retaining a bar. The exact 



significance of this mechanism is thus unclear but its 
effectiveness here is suggestive. In any case, it seems 
unlikely that it could produce a large bulge like that 
of, say, NGC 5866 (Sandage & Bedke 1994). 

Quinn et al. (1993) noted that subsequent ac- 
cretions by already thickened disks produced much 
smaller changes in thickness and dispersion than the 
initial accretion, so perhaps a bulge could be built 
up piecewise without damaging the disk excessively. 
Presumably this process would thicken any thin disk 
that had reformed after the previous encounter and 
the gas would gradually be depleted as some portion 
of it was driven to the centre each time. This leads to 
a picture in which a large disk galaxy, if it suffered sev- 
eral mergers with companions of non-negligible mass, 
might represent a nightmare scenario for anyone at- 
tempting to sort out its evolutionary history. It would 
have a thick disk composed of several populations of 
disk stars and stripped satellite stars. The disk would 
be relatively gas-poor with little present-day star for- 
mation and relatively hot with little present-day spi- 
ral structure. It would also likely have a dense nu- 
cleus composed of populations from several satellites 
and star formation episodes. SO galaxies are defined 
by their lack of spiral structure and they typically 
also have thicker disks, larger bulges, and less star 
formation than other disk galaxies (Sandage 1961; 
Sandage & Bedke 1994). Sa galaxies are defined by 
their tightly wrapped spiral structure, like that visi- 
ble in the last frame of Figure ^|. Thus mergers seem 
to help in all respects to push a spiral galaxy along 
the Hubble sequence toward Sa/SO, a contribution 
to galaxy evolution analogous to, but suitably scaled 
down from, the process of two spiral galaxies merging 
to form an elliptical. 

4.3. Is there a constraint on cosmology? 

Unfortunately, our results do not say much about 
the cosmological implications of companion accretion 
beyond the fact that disks are rather fragile. Because 
our satellites start at the edge of the disk (to save 
CPU time and beat the growth of noise), we do not 
really address sinking times from a cosmological per- 
spective. We find very short sinking times but a satel- 
lite at some distance out in the halo would slowly orbit 
many times before reaching the disk (Fig. |l3|). That 
is to say, the satellite is "accreted" by the halo long 
before it is accreted by the disk, which tells us that 
accretion rates derived from cosmology need to be in- 
terpreted carefully if they are to be constrained by 
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disk structure arguments. 

Figure ^ gives sinking times as a function of satel- 
lite mass and initial orbital radius for circular orbits 
decaying by dynamical friction in an isothermal halo 
(Binney & Tremaine 1987). At a glance one can see 
which kinds of objects will hit the disk within a Hub- 
ble time or, more importantly, within the age of the 
disk. Note that, while the full ranges of masses and 
orbital radii arc realistic for satellite galaxies, only a 
small region of the diagram has any relevance: not 
just any satellite can participate in the processes dis- 
cussed in this paper. Distant, low mass satellites take 
longer than the age of the disk to spiral in. The cal- 
culation is not applicable when the satellite is in the 
disk and any satellite there now must have started 
in the shaded region. Satellites whose mass is com- 
parable to the disk's mass do not qualify as "mi- 
nor" mergers. Evidently there is a roughly triangu- 
lar region in i? m j^ versus M sa ^ which we might call 
the "Zone of Interesting Parameters", where satel- 
lites must originate to be important. (There must 
also be a density limit below which a satellite will be 
disrupted before reaching the disk.) For illustration, 
the zone is bounded here by a 10 Gyr disk age, a 
20 kpc disk radius, and a 0.25^/^^. mass threshold 
for severe destruction, but naturally the zone is much 
more nebulous. In any case, the sinking times are 
merely rough estimates which will be improved only 
with better knowledge of the outer regions of galactic 
halos. Good galaxy-formation/cosmology simulations 
would be useful in this regard and would also help 
us determine the expected space, energy, and angu- 
lar momentum distributions of satellites around large 
galaxies. 

Specific comparison with the work of Toth & Os- 
triker (1992) is very difficult. In a vague way, the 
fact that the damage done to our disk is significant 
strengthens their result. However, their results are 
averaged over all possible orbital inclinations, making 
it impossible to directly relate our numbers to theirs. 
This averaging does not account for the fact that pro- 
grade orbits are more likely to merge within a Hubble 
time, nor does it account for the tendency of these or- 
bits to settle to low inclinations, depositing most of 
the vertical energy at large radii and favouring i ~ 
at smaller radii. Another difficulty is that their anal- 
ysis assumes the satellite's energy is deposited locally 
and its orbit decays slowly. These assumptions might 
be reasonably accurate for retrograde orbits but Fig- 
ures and 0, which illustrate the global response of 



the disk and the rapid orbital decay respectively, show 
that our prograde mergers violate them. 

Toth & Ostriker set a limit on accretion by start- 
ing with an absolutely cold disk and heating it with 
giant molecular clouds, which gives Q — 0.90 and 
a scaleheight h = 218 pc. They then attribute the 
difference between these values and the present thin 
disk values to accretion. The assumptions that total 
heating and thickening scale with total accreted mass 
and that they are the same regardless of the initial 
heat and thickness are not strictly correct. Quinn 
et al. (1993) found that a second satellite impact on 
an already heated disk had a much less dramatic ef- 
fect than the first: two satellites of O.Q5Mdisk do not 
equal one satellite of O.lMdisk- In fact, in the limit 
where the satellite merger is finely subdivided into 
smaller events, this goes over to an adiabatic pro- 
cess of slow, continuous accretion or is shut off alto- 
gether if the objects are too distant to sink within 
a Hubble time (Fig. |l3|). Also the assumption of a 
low pre-merger Q means that the disk is unstable to 
global structure formation and will experience addi- 
tional "initial" (non-merger) heating. Finally, when 
Toth & Ostriker take account of mass loss by the 
satellite (which they otherwise treat as rigid), they 
use the density ratio approximation, which we have 
found to be inadequate (§ |3.3| ). Some of these as- 
sumptions and approximations actually lead Toth & 
Ostriker to underestimate their damage; the implica- 
tions of others are less clear. Overall, it is not clear 
that they have succeeded in placing a limit on ac- 
cretion, nor that our results strengthen any of their 
conclusions. 

The bottom line is that two things are needed for a 
limit on Clo: a limit on accretion and a connection be- 
tween accretion and cosmology. It is not certain that 
we have either, yet. Meanwhile, it is worth noting 
that if gas can reform a thin disk after an encounter, 
then the presence of a thin disk does not necessarily 
mean there have been no encounters. 

5. Conclusions 

Our fully self-consistent simulations of the accre- 
tion of a companion galaxy by a large disk galaxy 
allow us to examine the full interaction between the 
satellite, disk, and dark matter halo. We find that 
500000 particles is sufficient to eliminate numerical 
relaxation over the timescale of the prograde merger, 
but mergers not strongly coupled to the disk (retro- 
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grade, polar) still face competition from noise. 

That the disk should dominate the interaction 
could have been predicted from the coupling between 
the satellite's orbit and the orbits of disk stars. The 
strong global response of the disk depends on a satel- 
lite orbit which is both prograde and circular. Either 
high inclination or high eccentricity would curtail the 
resonance considerably. It remains to be seen whether 
the "typical" satellite orbit becomes circular by the 
time it reaches the optical disk, but our initial condi- 
tions ought to be viewed as a somewhat special case. 
Future stellar dynamical studies of satellite accretion 
should examine retrograde and eccentric orbits. We 
expect them to find that the halo contributes a larger 
fraction of the torque and that less of the satellite 
survives to the galactic centre. 

With its high density and the rapid decay of its pro- 
grade orbit, our satellite easily survives the encounter. 
Half of its mass arrives in the central few kiloparsecs 
of the disk, causing significant heating there. After 
the merger, the satellite core is a small, bulge-like 
entity in the galactic centre while the galactic disk is 
rather like the thick disk component of the Milky Way. 
Our particular encounter produced only about 60% 
thickening at the solar circle but enough heating so 
that the final (ctr, <j z )q ~ (48,42,38) km/s, very 
close to the observed values for o$ and er z . Our <tr 
and asymmetric drift are lower than those observed 
but a more eccentric satellite orbit might boost them 
a bit. Eventually, these effects all leave the galaxy 
looking like an earlier Hubble type (Sa in this case, 
since there is spiral structure) , supporting the picture 
of merger-induced galaxy evolution. 
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Fig. 1. — Initial radial disk structure. Particles were 
evenly binned into 100 cylindrical shells and quanti- 
ties are thus averaged over azimuth. Clockwise from 
the upper left, the panels display surface density, 
disk thickness, rotation speed and asymmetric drift, 
and velocity ellipsoid. Disk thickness is represented 
by (| z |) which is related to zq (where zo is approxi- 
mately twice the exponential scaleheight; see eq. [I]) by 
= zo In 2, but only to the extent that this trans- 
formation is applicable. Scaleheights are only mean- 
ingful insofar as the particle distribution fits the rele- 
vant functional form. is unambiguous because no 
quality-of-fit indicator is needed for its interpretation. 



Fig. 2. — Early equilibrium disk structure. This fig- 
ure shows for an isolated galaxy the same quantities 
as Figure |l| but at time t — 63 Myr, after the galaxy 
has come into equilibrium. For comparison, the t = 
structure is superposed with dotted curves. Note in 
particular the 20% drop in the central velocity dis- 
persions. 



Fig. 3. — Structure evolution at the solar radius in 
isolated galaxy models. These are plotted on the same 
vertical scales as the corresponding panels in Figures [l] 
and [2] to facilitate visual comparison. The solid curves 
represent the 450 000-particle simulation, dotted rep- 
resent N = 225 000, dot-dashed N = 90 000, and 
dashed N — 45 000. This figure shows that we have 
defeated numerical relaxation noise over the gigayear 
timcscale required for our prograde mergers. Note 
that there is an initial transient because the disk does 
not start in perfect equilibrium (as shown in Fig. |^ 
and discussed in § 2.1). Once they have settled, the 
large runs remain essentially unchanged until their 
dispersions begin to creep up at t ~ 3 Gyr. 
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Fig. 4. — Growth of m = 2 (bar) mode in isolated 
galaxy models. This figure displays the growth of 
bars in the same simulations shown in Figure || and 
also in a run with 45 000 disk particles and a rigid 
halo. The bars set in at later times in larger simula- 
tions because they are seeded by random clumping, 
the amplitude of which decreases when more particles 
are used. That the halo is the source of these seeds 
in the self-consistent simulations is made clear by the 
lack of growth in the rigid halo case. The largest 
run suppresses the bar long enough for our prograde 
mergers to complete, but not our retrograde mergers. 
This illustrates the necessity of using very large N, 
particularly in the halo, to smooth out the potential 
in self-consistent simulations. 

Fig. 5. — Face-on and edge-on views of the disk and 
satellite particles at equal intervals of ~ 125 Myr, 
starting at t = 0. The disk's global response to the 
satellite is quite apparent in the face-on panels while 
the thickening and warping of the disk are apparent 
in the edge-on view. The global tilt is removed before 
further analysis by a rotation which aligns the total 
angular momentum vector of the disk particles with 
the z-axis. The satellite core arrives at the centre in 
the penultimate frame. 

Fig. 6. — Satellite orbit. Clockwise from the up- 
per left, the panels show the cylindrical radius ver- 
sus time, y versus x (corresponding to the "face-on" 
view in Figure ||), altitude versus radius, and alti- 
tude versus time. Our satellite completes fewer than 
two orbits before intersecting the solar circle and our 
merger is all over by t = 1 Gyr. Note in particular 
that the satellite settles into a low-inclination orbit 
while it is still at a large radius. (For comparison, the 
dotted lines in the lower right panel show the initial 
inclination, i = 30°.) 



Fig. 7. — The upper panel shows the torques acting 
on the satellite core due to the disk and halo (in sim- 
ulation units). The disk is responsible for about 75% 
of the total torque integrated over the duration of the 
merger. The lower panel superposes the disk torque 
on the bar mode and orbital decay curves to make 
their close relationship clear. 



Fig. 8. — The upper panel compares the decay of the 
satellite orbit in four simulations of the same orbit 
(our main, prograde encounter) covering an order of 
magnitude in size to test convergence. The consis- 
tency is very encouraging. Also shown is the decay of 
the same orbit when a rigid halo is used. The qual- 
itative behaviour is the same but, because the satel- 
lite has nothing to interact with when it is far from 
the disk plane, the interval preceding the knee in the 
curve is longer. Once the satellite has settled into the 
disk plane, the decay rates in the rapid sinking phase 
are about the same because the disk dominates. This 
is also true of the Quinn et al. (1993) simulations 
except that the truncation of their disk makes the in- 
terval before rapid sinking even longer, about 2 Gyr. 
The lower panel compares the decay of our fiducial 
orbit with that of other orbits: retrograde with 30° 
inclination, polar, and coplanar. Clearly, the polar 
case is not intermediate between the prograde and 
retrograde cases but is essentially retrograde. This 
demonstrates the significance of the strong disk cou- 
pling for prograde orbits: it leads to a factor of three 
reduction in the sinking time (note the different scal- 
ing of the time axes in the upper and lower panels). 

Fig. 9. — Post-merger disk structure. The solid 
curves show the disk in our encounter and the dashed 
curves show the disk in the isolated galaxy. (This 
comparison accomodates both the initial transient 
and any common numerical relaxation effects.) Both 
are shown at t = 1.19 Gyr, which leaves time for 
things to settle down, although not much changes be- 
tween t = 1.0 Gyr and t = 1.2 Gyr. After this time, 
the only significant changes are associated with the 
bar's vertical instability (see § ||). Note the greater 
radial extent of the disk which underwent the merger 
(solid curves), indicating conversion of satellite orbital 
energy to disk potential energy. 

Fig. 10. — Disk structure evolution at the solar cir- 
cle (i?0 = 8.0 kpc). The top panel illustrates the 
thickening of the disk while the centre panel shows 
the velocity ellipsoid and asymmetric drift. The lower 
panel repeats the orbital decay curve from Figure ^ so 
that features in the disk structure can be temporally 
matched up with the satellite's location. In partic- 
ular, the dotted, vertical line shows the time when 
the satellite crosses the solar circle. Note that most 
quantities rise abruptly when the satellite crosses Rq . 
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Fig. 11. — Integrated mass distribution and surface 
density of the satellite remnant. That almost half the 
satellite escapes tidal stripping and forms a compact 
central element is apparent in both panels. The up- 
per panel also shows the retrograde case (in which 
the satellite is nearly disrupted) for comparison. The 
time is 1.19 Gyr. In the lower panel, the satellite's 
surface density curve reveals that the core has not 
quite finished sloshing around the centre of the disk. 
The curve for the disk and satellite combined shows 
that the satellite core is significant enough to notice- 
ably enhance the central brightness of the disk. 

Fig. 12. — Structure of disk + satellite combined. The 
solid curves show the mass- weighted average structure 
of the luminous material (disk+satellite). For com- 
parison, the disk is shown with dashed curves and 
the satellite remnant with dotted curves. The time is 
1.19 Gyr. Only in the inner two kiloparsecs is there 
a sufficient density of satellite material to cause the 
combined structure to deviate substantially from the 
disk structure. Outside the core, the satellite rem- 
nant is like a hot, thick, flared disk, although most of 
its "thickness" at large radii (R > 15 kpc) actually 
represents tilting because that material is still in the 
original, inclined orbital plane. 

Fig. 13. — Sinking isochrones for an isothermal halo 
(based on eq. 7-27 in Binney & Tremaine (1987)). 
This figure provides estimates of the sinking time for 
a satellite of given mass starting in a circular orbit 
of given radius. The only objects which are signifi- 
cant for the present study are those which originate in 
the shaded, triangular region and are dense enough to 
survive disruption and reach the disk. All other satel- 
lites either sink too slowly or are already in the disk or 
are so massive they will destroy the disk. There may 
also be an upper cutoff if the tidal radius of a typi- 
cal halo is much smaller than 100 kpc. For reference, 
our satellites all originate at the position marked with 
the cross, just above the lower boundary (for compu- 
tational expedience). That the diagram seems to pre- 
dict the right sinking time for the prograde encounter 
is just an accident. All our simulations started at 
the cross but their sinking times spanned 0.6 Gyr to 
3+ Gyr; the uncertainties in the isochrone positions 
are also quite large. 



